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Abstract 

We adopt a factorized QCD formalism to describe the transverse momentum distribution of low- 
mass lepton pairs produced in pp collisions, when the pair transverse momentum Qj » Q, 
with the pair's invariant mass Q as low as Q ~ Aqcd- We extend this formalism to dilepton 
production in AA collisions by including the nuclear-dependent power correction due to parton 
multiple scattering. 



1. Introduction 

Dilepton production in ultrarelativistic heavy-ion collisions is among the richest sources of 
information on the features of the produced strongly-interacting medium. Since leptons do not 
participate in the strong interaction directly, dilepton production could probe the entire space- 
time evolution of the expanding hot dense system without the interference of final-state interac- 
tions. An enhanced dilepton yield from the thermal emission of the hot medium could be a signal 
of the deconfining phase transition for the strongly interacting matter fl]] . 

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) has recently mea- 
sured the transverse momentum distribution of electron-positron pairs, with the pairs' invariant 
mass as low as 100 < Q < 300 MeV, in both pp and AA collisions. For Qj < 3 GeV they 
discovered a strong excess of dileptons in AA collisions in comparison with the binary-collision 
scaled pp result ||2|]. 

To have a reliable prediction for the dilepton yield in this low-mass regime Qj » Q ~ Aqcd 
is a challenge for perturb ative QCD, even in the case of pp collisions. We argue that the distribu- 
tion at large Qj can be systematically factorized into universal parton-to-dilepton fragmentation 
functions, parton distributions, and perturbatively calculable partonic hard parts evaluated at a 
short distance scale ~ 0(1 / Qt) 13]- We introduce a model for the input dilepton fragmentation 
functions at a scale hq ~ 1 GeV, which are then evolved perturbatively to scales relevant at RHIC. 
Using the evolved fragmentation functions, we evaluate the transverse momentum distributions 
in pp collisions and are able to describe the PHENIX data in pp collisions quite well. 

We further extend our formalism to dilepton production in AA collisions by considering both 
the leading-power nuclear effects from nuclear parton distribution functions (nPDFs), and the 
nuclear-dependent power corrections generated by parton multiple scattering. We find that the 
shadowing of nPDFs suppresses the dilepton yield while the nuclear power correction enhances 
it. But, the nuclear power correction alone is not sufficient to explain the large dilepton excess 
observed by PHENIX. We also explore the possibility that in AA collisions the hot dense medium 
is produced before the hard collision takes place and discuss the impact this would have on 
generating the dilepton excess. 

Preprint submitted to Nuclear Physics A September 23, 2009 



2. QCD resummation formalism for low-mass dilepton yield in pp collisions 



The cross section for Drell-Yan type inclusive dilepton production in hadronic collisions, 
A(Pa)+B(Pb) — > > ^£~(0)+X, can be expressed in terms of the cross section for producing 
a virtual photon that decays into the observed dilepton: 
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When both physically measured momentum scales Q and Qj are much larger than Aqcd, the 
cross section for producing the virtual photon can be factored systematically in QCD perturbation 
theory as J3l 
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where Yia)> mns over au> parton flavors and d& F a "^,^ x l dQ T dy represents corresponding per- 
turbatively calculable short-distance hard parts. When Qj » Q, high order hard parts could 
receive large corrections in powers of a s \n{Q^ / Q 2 ) caused by the radiation of partons along the 
direction of the low-mass dilepton. Such large logarithmic corrections can be systematically re- 
summed into the parton-to- virtual photon fragmentation functions, Df^y> fl. The perturbative 



hard parts, d&^^,^ x / dQ T dy, can be re-organized into two terms as |5| 
c/o- Pert dfr Dk dfr Fms 
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where the superscripts "Dir" and "Frag" represent the "direct" and the "fragmentation" contribu- 
tion, respectively. The fragmentation contribution can be further factorized as |5[], 
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with all fragmentation logarithms resummed into Df-, y > and the perturbative dcr a i,^ c x/dp 2 T dy 
evaluated at the distance scale 1/Qt- The "Dir" term is a perturbative difference between the 
"Pert" and "Frag" terms when both are expanded in powers of a s . 

We have argued in Ref. [3] that the factorization formalism for the dilepton yield at high 
Qj is also valid for lepton pairs with invariant mass Q ~ Aqcd- However, in this case the 
fragmentation functions, D c ^ y >, are no longer fully perturbative but introduce non-perturbative 
hadronic contributions to the dilepton yield. We proposed to account for these by adding a 
hadronic component to the fragmentation functions at a factorization scale /jq ~ 1 GeV 0]: 



D^fcpg; g 2 ) = Df^\z,p.l; Q 2 ) + kD^hD = 
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where D^ D » represents the leading order QED component and the size of the hadronic compo- 
nent is fixed by a fitting parameter k. Following we further assume that V = p, fp/4-jr = 2.2, 
and Df^ p w D/-^. In the top left figure of Fig. [T] we compare our calculation to the PHENIX 
dilepton data in pp collisions [2] under two assumptions: k — 1 (solid line) and k — (dotted 
line). Our calculation based on the factorization formalism describes the data rather well. The 
hadronic component is very important when Qj < 3 GeV. 
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3. Nuclear modification in AA collisions 



There could be many interesting sources of nuclear modification for dilepton production in 
AA collisions. For the leading-power factorization formalism, one could have a universal nuclear 
effect from the nPDFs. We can study this effect by replacing the proton PDFs in Eq. (|2|l by the 
nPDFs, ff -> [Z ■ fj ,/A + (A - Z) ■ f" /A ] /A. Using the EPS08 nPDFs H, we plot our calculated 
cross sections with k — 1 (dashed curves) in Fig. Q] and compare to the PHENIX data in Au+Au 
collisions yfl. For comparison, we also show the binary-collision scaled pp results as dotted 
curves. The cross section with nPDFs is larger than the scaled pp results when Qj > 3 GeV 
due to the antishadowing effect and is consistent with the data. However, the shadowing effect 
in nPDFs leads to a suppression of the dilepton yield when Qj < 3 GeV, where the data shows 
a very strong enhancement over the scaled pp results. 
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Figure 1: Invariant cross section for low-mass dilepton production with 100 < Q < 300 MeV as a function of the pair's 
transverse momentum at ySjviv = 200 GeV and y = 0. Top left: in pp collisions with k = 1 (solid) and k = (dotted). 
Remaining three plots: in Au+Au collisions with different centralities. The dashed and solid lines are evaluated with 
nPDFs, and nPDFs plus power correction, respectively. The dotted lines are for binary-collision scaled pp cross section. 
The data are from Ref. Ql- 

In nuclear collisions, the high-2r lepton pairs could also be produced by multiple scatter- 
ing. When Q and Qj are small, such initial-state multiple scattering could lead to medium-size 
enhanced power corrections to the cross section and significantly increase the yield of low-mass 
lepton pairs [8]. The leading power correction to the dilepton yield is proportional to the twist-4 
quark-gluon and four-gluon correlation functions, T q /A(x) and T g /A(x) Q8[] . We evaluate the power 
correction by adopting the usual model T qg /A(x) = A 2 A 4 ^f q / A (x), with f q l g (x) denoting as be- 
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fore the proton PDFs in nuclei and A 2 w 0.01 GeV 2 fixed by data [9]. We plot our predictions for 
the dilepton yield resulting from inclusion of this power correction as solid lines in Fig. [Q We 
find that the power correction could enhance the yield by 70 - 90% at low Qj. However, the cal- 
culated cross section including the power correction is still significantly below the experimental 
data at low Qj. 

Our calculation indicates that the power correction from initial-state parton multiple scatter- 
ing in cold nuclear matter cannot explain the dilepton excess observed by PHENIX. As argued 
m Ref. [10], the initial density of the hot medium produced in central Au+Au collisions at RHIC 
could be as large as 30 times that of cold nuclear matter. Should this hot medium (or a part of it) 
be formed before the hard collision that produces the high-Q^ lepton pair, the rate for multiple 
scattering, and thus the dilepton yield from multiple scattering, could be significantly enhanced 
by the higher medium density. To fit the PHENIX data, we would need a hot medium density 
that is an order of magnitude higher than the cold nuclear density. 

4. Conclusion 

We have investigated the production cross section for lepton pairs in the regime Qj » Q ~ 
Aqcd hi both pp and AA collisions. We have argued that the cross section for \\\g\\-Qj low-mass 
dileptons can be treated in perturbative QCD in the same way as the cross section for direct 
photon production. Our calculation is consistent with the RHIC data for pp collisions. We have 
also calculated the yield of low-mass lepton pairs in AA collisions by including nuclear effects 
from nPDFs and parton multiple scattering in cold nuclear matter. We have found that multiple 
scattering enhances the dilepton yield at low Qj, while the nPDFs suppress it. Although the net 
effect enhances the dilepton yield, the enhancement is far too small to explain the dilepton excess 
observed by the PHENIX collaboration. The data requires us to consider additional sources of 
low-mass lepton pairs in AA collisions. The thermal emission of low mass lepton pairs from 
the hot medium could be an important source U]. However, because of the high transverse 
momentum, the thermal emission alone might not be able to generate the additional order of 
magnitude more lepton pairs at Qj > 1 GeV needed to explain the data. 
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